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(Dated: January 9, 2007) 

From an analysis of the decay B° — > J/ip(j) we obtain the width difference between the light and 
heavy mass eigenstates, Ar = (T L - F H ) = 0.17 ± 0.09 (stat) ±0.02 (syst) ps" 1 and the CP- 
violating phase 4> s — —0.79 ± 0.56 (stat) l"o'oi ( s y s *0- Under the hypothesis of no CP violation 
(4> s = 0), we obtain 1/T = t(B°) =1.52 ±0.08+°;°,3 ps and Ar = 0.12+q;J 8 ± 0.02 ps" 1 . The data 
sample corresponds to an integrated luminosity of 1.1 fb~ 4 accumulated with the DO detector at 
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the Tevatron. 

PACS numbers: 13.25.Hw, 11.30.Er 

In the standard model (SM), the light (L) and heavy 
(H) eigenstatcs of the mixed B® system are expected 
to have a sizeable mass and decay width difference, 
AM = M H - Ml and Ar = T L - T H . The CP-violating 
phase, defined as the the relative phase of the off-diagonal 
elements of the mass and decay matrices in the B® - B s 
basis, is predicted to be small. Thus, to a good approxi- 
mation the two mass eigenstates are expected to be CP 
eigenstates. New phenomena may alter the CP-violating 
mixing phase </> s , leading to a reduction of the observed 
Ar compared to the SM prediction Q Ar SA /: Ar = 
Ar^M x cos 4> s - While the mass difference has recently 
been measured to high precision 0, Q > the CP-violating 
phase remains unknown. 

The decay B® — > J / ipcp, proceeding through the quark 
process b — > ccs, gives rise to both CP-even and CP- 
odd final states. It is possible to separate the two CP 
components of the decay B® — ► J/ip<fi, and thus to mea- 
sure the lifetime difference, through a study of the time- 
dependent angular distribution of the decay products of 
the J/ip and <fi mesons. Moreover, with a sizeable life- 
time difference, there is sensitivity to the mixing phase 
through the interference terms between the CP-even and 
CP-odd waves. 

In Ref. [3| we presented an analysis of the decay chain 
B° -> J/ip(j>, J/tp -> -» K + K~ based on the 

first w450 pb" 1 of pp data at a center-of-mas energy of 
1.96 TeV collected with the DO detector In that anal- 
ysis, we extracted three parameters characterizing the B® 
system and its decay B° B — ► J/ipcj>: the average lifetime, 
t = 1/r, where r ee (T H + T L )/2; AT/T; and the relative 
rate of the decay to the CP-odd states at time zero. Here 
we present new results, based on a two-fold increase in 
statistics. In addition to r and Ar, we extract for the 
first time the CP-violating phase <p s - We also measure 
the magnitudes of the decay amplitudes, and their rela- 
tive phases. 

The data, collected between June 2002 and January 
2006, correspond to an integrated luminosity of 1.1 fb _1 . 
The selected events include two reconstructed muons of 
opposite charge, with a transverse momentum greater 
than 1.5 GeV and pseudorapidity \r/\ < 2. Each muon 
is required to be detected as a track segment in at least 
one of the three layers of the muon system, and to be 
matched to a central track. One muon is required to have 
segments both inside and outside the toroid magnet. We 
require the events to satisfy a muon trigger that does not 
include a cut on the impact parameter. 

To select the B® candidate sample, we set the mini- 
mum values of momenta in the transverse plane for B® , 
</>, and K meson candidates at 6.0 GeV, 1.5 GeV, and 
0.7 GeV, respectively. J/ip candidates are accepted if 
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FIG. 1: The invariant mass distribution of the (J/ip,<f)) sys- 
tem for candidates. The curves are projections of the 
maximum likelihood fit (see text). 

the invariant mass of the muon pair is in the range 2.9 
- 3.3 GeV. For events in the central rapidity region (an 
event is considered to be central if the higher px muon 
has 1 77^1 1 < 1), we require the transverse momentum of 
the J/ip meson to exceed 4 GeV. Successful candidates 
are constrained to the world average mass of the J/ip me- 
son Q . Decay products of the <j> candidates are required 
to satisfy a fit to a common vertex and to have an invari- 
ant mass in the range 1.01 - 1.03 GeV. We require the 
{J/ip, 4>) pair to be consistent with coming from a com- 
mon vertex, and to have an invariant mass in the range 
5.0 - 5.8 GeV. In the case of multiple <p meson candidates, 
we select the one with the highest transverse momentum. 
Monte Carlo (MC) studies show that the pr spectrum of 
the 4> mesons coming from B® decay is harder than the 
spectrum of a pair of random tracks from hadronization. 
We define the signed decay length of a B° s meson L^ y as 
the vector pointing from the primary vertex to the de- 
cay vertex projected on the B® transverse momentum. 
To reconstruct the primary vertex, we select tracks with 
Pt > 0.3 GeV that are not used as decay products of 
the B® candidate, and apply a constraint to the aver- 
age beam spot position. The proper decay length, ct, is 
defined by the relation ct = L^ y ■ Mb»/pt where Mgo 
is the measured mass of the B® candidate. The distri- 
bution of the proper decay length uncertainty a(ct) of 
Bg mesons peaks around 25 /zm. We accept events with 
a(ct) < 60 fim. The invariant mass distribution of the 
accepted 23343 candidates is shown in Fig.[TJ The curves 
are projections of the maximum likelihood fit, described 
below. The fit assigns 1039±45 (stat) events to the B® 
decay. 
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We perform a simultaneous unbinned maximum likeli- 
hood fit to the proper decay length, three decay angles, 
and mass. The likelihood function £ is given by: 



N 

£ = JJ[/siff-^sig + (1 — /siff)-^6cfc]> 



(1) 



where TV is the total number of events, and f S i g is the 
fraction of signal in the sample. The function T\ ig de- 
scribes the distribution of the signal in mass, proper de- 
cay length, and the decay angles, and Tl ck is the prod- 
uct of the background mass, proper decay length, and 
angular probability density functions. Background is di- 



J 



vided into two categories. "Prompt" background is due 
to directly produced J/tp mesons accompanied by ran- 
dom tracks arising from hadronization. This background 
is distinguished from "non-prompt" background, where 
the J/ip meson is a product of a B hadron decay while 
the tracks forming the <j> candidate emanate from a multi- 
body decay of the same B hadron or from hadronization. 

The time evolution of the angular distribution of the 
products of the decay of flavor untagged B° mesons, i.e., 

summed over B® and B s , expressed in terms of the linear 
polarization amplitudes A x and their relative phases Si 
is0: 



d 3 r(t) 



dcosO dip (icos-0 



oc 2\A (0)\ 2 T+ cos^(l-shr 6>cosV)+sin 2 ^{\A {1 (0)| 2 T+ (1-sin 2 6>sin>) + |,4 ± (0)| 2 T_ sin 2 9} 



+ ^sin2^|Ao(0)||A||(0)|cos(J 2 -Si) T+ sin 2 6> sin 2cp 
v2 



f { —^=1^0(0)11^(0)1 cos5 2 sin 2tp sin 29 cos if 



-\A\\(0)\\Ax (0) | cos 6i sin 2 tp sin 29 sin i {e~ TH 1 - e~ TL *) sin S 



(2) 



r 



-r H t\ 



where T + /_ = \ ((1 ± cos0 s )e Vl1 + (1 T cos^ s )e 

In the coordinate system of the J/ ip rest frame (where 
the 4> meson moves in the x direction, the z axis is 
perpendicular to the decay plane of <f> — > K + K~ , and 
p y (K + ) > 0), the transversity polar and azimuthal an- 
gles (9, if) describe the direction of the and ip is the 
angle between p(K + ) and —p(J/ip) in the <f> rest frame. 

We model the acceptance in the three angles by fits 
using polynomial functions, with parameters determined 
using Monte Carlo simulations. We have used the 
SVV.HELAMP model in the EvtGen generator 0], in- 
terfaced to the Pythia program [§]. Simulated events 
were reweighted to match the kinematic distributions ob- 
served in the data. 

The lifetime distribution shape of the background is 
described as a sum of a prompt component, simulated 
as a Gaussian function centered at zero, and a non- 
prompt component, simulated as a superposition of one 
exponential for the negative ct region and two exponen- 
tials for the positive ct region, with free slopes and nor- 
malization. The mass distributions of the backgrounds 
are parametrized by first-order polynomials. The dis- 
tributions in the transversity polar and azimuthal an- 



gles are parametrized as 1 + X? x cos 2 9 + X^ x cos 4 9 and 
1 + cos(2yj) +Y2x cos 2 (2</?), respectively. For the back- 
ground dependence on the angle -0, we use the function 
1 + Z 2x cos 2 (-0). We also allow for a background term 
analogous to the interference term of the CP-even waves, 
with one free coefficient. For each of the above back- 
ground functions we use two separate sets of parameters 
for the prompt and non- prompt components. 

Our results for the hypothesis of CP conservation and 
for the case of free <p s , are presented in Table HI For the 
CP-violating phase, which has a four-fold ambuguity dis- 
cussed below, the fit value closest to the SM prediction 
of -0.03 Q is <j> s = -0.79 ± 0.56. Figures - show 
the fit projections on the angular distributions and the 
proper decay length. Figure [6] shows the Aln(£) = 0.5 
error ellipse contour (corresponding to the confidence 
level of 39%) in the plane (Ar, 4> s ). As seen from Eq. [2J 
the sign of sin (f> s is reversed with the simultaneous re- 
versal of the signs of cos (Si and cos<5 2 . For the case 
cos<5i < and cos^ > 0, our measurement correlates 
two possible solutions for <p s with the two signs of Ar: 
4> s = -0.79±0.56, Ar > 0, and 4> s = 2.35±0.56, Ar < 0. 
For the case cos Si > and cos 82 < the two solutions 
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are (p a - 
Ar < 0. 



0.79±0.56, Ar > 0, and 



-2.35±0.56, 



TABLE I: Maximum likelihood fit results. Sign ambiguities 
are discussed in the text. 



Observable 


CP conserved 


free <j> s 


Ar (ps- 1 ) 


12+ 008 
u - lz -o.io 


n 1 7 +o.09 

""■■'-0.09 


b = r 


r 9 +0.08 


1.49 ±0.08 




= 


-0.79±0.56 


L4 () (0)! 2 -|^||(0)| 2 


0.38±0.05 


0.37±0.06 


A ± (0) 


0.45±0.05 


0.46±0.06 


Si — s 2 


2.6±0.4 


2.6±0.4 


8i 




3.3±1.0 


52 




0.7±1.1 
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• Data 

— Total Fit 
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— Background 




PIG. 3: The transversity asimuthal angle distribution for 
the signal-enhanced subsample: ct/a(ct) > 5 and signal mass 
range. The curves show: the signal contribution, dotted (red) ; 
the background, light solid (green); and total, solid (blue) 
[color online]. 
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PIG. 2: The transversity polar angle distribution for the 
signal-enhanced subsample: ct/cr(ct) > 5 and signal mass 
range. The curves show: the signal contribution, dotted (red); 
the background, light solid (green); and total, solid (blue) 
[color online]. 
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PIG. 4: The i/j angle distribution for the signal-enhanced 
subsample: ct/a(ct) > 5 and signal mass range. The curves 
show: the signal contribution, dotted (red); the background, 
light solid (green); and total, solid (blue) [color online]. 



We perform a test using pseudo-experiments with sim- 
ilar statistical sensitivity, generated with the same pa- 
rameters as obtained in this analysis under the condition 
of no CP violation. When fits allowing for CP violation 
are performed, w 50% of the experiments have a fitted 
cos(c/) s ) less than the measured value. About 80% of ex- 
periments have the statistical uncertainty of cf> s greater 
than that for data. 

We verify the procedure by performing fits on MC 
samples passed through the full chain of detector sim- 
ulation, event reconstruction, and maximum likelihood 
fitting. We assign systematic uncertainties due to the 
statistical precision of this procedure test. We also re- 
peat the fits to the data with the parameters describing 
the acceptance varied by ±ltr. 



Uncertainties from the data processing reflect the sta- 
bility of the results with respect to different versions of 
the track and vertex reconstruction algorithms. The "in- 
terference" term in the background model accounts for 
the collective effect of various physics processes. How- 
ever, its presence may be partially due to the detector 
acceptance effects. Therefore, we interpret the difference 
between fits with and without this term as a systematic 
uncertainty associated with the background model. Ef- 
fects of the imperfect detector alignment are estimated 
using a modified geometry of the the silicon microstrip 
tracker, with silicon sensors moved within the known 
uncertainty. The effects of systematic uncertainties are 
listed in Table HH 

From a fit to the CP-conserving time-dependent angu- 
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1.1 fb' B° s -> J/\|/<|) • Data 

i Mass 5.26 - 5.46 GeV Total Fit 

■ Total Signal 
CP-even 
CP-odd 
— Background 



0.25 0.3 
Ct (cm) 

FIG. 5: The proper decay length, ct, of the B® candidates 
in the signal mass region. The curves show: the signal con- 
tribution, dashed (red); the CP-even (dotted) and CP-odd 
(dashed-dotted) contributions of the signal, the background, 
light solid(green); and total, solid (blue) [color online]. 



TABLE II: Sources of systematic uncertainty in the results of 
the analysis of the decay B® — > J/ip(j>. 



_0.5 
in 

5-0.4 

1— 

u0.3r> 1 

< 

0.2E 

0.1 
-0 
-0.1 
-0.2 
-0.3 
-0.4 



-°- 5 -5 



D0 

B°- 



1 .1 fb 

• J/\)/<j) 



— SM 

Ar = Ar SM x|cos((yi 



-4 



-2 -1 







1 



3 4 5 
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FIG. 6: The Aln(£) = 0.5 contour (error elipse) in the plane 
(AT, <f>a) for the fit to the B® — > J/ip4> data. Also shown is 
the band representing the relation AF = ATsm X |(cos(^ s )|, 
with Ar SA / = 0.10 ± 0.03 ps" 1 [jjj- The 4-fold ambiguity is 
discussed in the text. 



lar distribution of the untagged decay B® — > J/ip(f), we 
obtain the average lifetime of the B® system, t(B®) = 
1.52 ±0.08 (stat) 1^03 ( s y st ) P s and thc width difference 
between the two mass eig enstates, AT = 0.12^°;°^ (stat) 
±0.02 (syst) ps" 1 . 



Allowing for CP violation in B® mixing, we provide the 
first direct constraint on the CP-violating phase, <j) s = 
-0.79 ±0.56 (stat) to.lt (syst). 



Source 


ct(B») 
/an 


Ar 

ps" 1 


Rx 


4>s 


Procedure test 


±2.0 


±0.02 


±0.01 




Acceptance 


±0.5 


±0.001 


±0.003 


±0.01 


Reco. algorithm 


-8.0,±1.3 


±0.001 


±0.01 


-0.01 


Background model 


±1.0 


±0.01 


-0.01 


±0.14 


Alignment 


±2.0 








Total 


-8.8, ±3.3 


±0.02 


±0.02 


-0.01, ±0.14 
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